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ABSTRACT 

We present Berkeley Illinois Maryland Association (BIMA) millimeter interferometer observations 
of giant molecular clouds (GMCs) along a spiral arm in M31. The observations consist of a survey 
using the compact configuration of the interferometer and follow-up, higher-resolution observations on 
a subset of the detections in the survey. The data are processed using an analysis algorithm designed 
to extract GMCs and correct their derived properties for observational biases thereby facilitating 
comparison with Milky Way data. The algorithm identifies 67 GMCs of which 19 have sufhcient signal- 
to-noise to accurately measure their properties. The GMCs in this portion of M31 are indistinguishable 
from those found in the Milky Way, having a similar size-line width relationship and distribution of 
virial parameters, confirming the results of previous, smaller studies. The velocity gradients and 
angular momenta of the GMCs are comparable to the values measured in M33 and the Milky Way; 
and, in all cases, are below expected values based on the local galactic shear. The studied region of 
M31 has a similar interstellar radiation field, metallicity, Toomre Q parameter, and midplane volume 
density as the inner Milky Way, so the similarity of CMC populations between the two systems is not 
surprising. 

Subject headings: Galaxies: individual (Andromeda) — galaxies: ISM — ISM: clouds — radio lines: 
ISM 



1. INTRODUCTION 

As the instrumentation for millimeter-wave telescopes 
improves, it becomes progressively more straightforward 
to study individual molecular clouds in other galaxies. 
Recent studies of Local Group galaxies have surveyed 
large numbe rs of molecular cloud s in the Large Magel- 
lanic Cl oud iMizuno et al.ll2001bD . the Small Magellanic 
Cloud l|Mizuno et all l2001a|) . M33 (jEngargiola et alJ 
|2003|) , and a bev y of Local Group dwarf galaxies 
(e.g. IWilsonl IT99I iTavlor et al.lll999D . These recent 
studies explore the nature of star formation on galac- 
tic scales by studying the properties of giant molecu- 
lar clouds (GMCs, M > 10^ Mq) throughout their host 
galaxies. Such GMCs contain the majority of the molec- 
ular mass in the Milky Way's ISM and are responsibl e 
for most of the star formation in the Galaxy l)Blitdll99^ . 

The Andromeda Galaxy (M31) is the second largest 
disk galaxy in the Local Group, after the Milky Way, 
and it subtends over 2 deg^ on the sky. Its proximity 
(770 kpc, Jrccdman \: .Madorc 1990) makes it an excel- 
lent target for studying extragalactic molecular clouds. 
Numerous surveys of CO emission have been conducted 
over a portion of M31 and a comprehensiv e list of the 24 
CO st udies published up to 1999 is given in lLoinard et al.l 
l)1999D . This extensive list of surveys can be supple- 
ment ed with a few major studie s that have occurred since 
then. iSheth et alJ l|2000l 120061) used the BIMA millime- 
ter interferometer to study a 3' field in the outer region 
of the galaxy {Rgai — 12 kpc) and find 6 molecular com- 
plexes similar to those found in the Milky Way. An ex- 
tensive survey covering the entirety of the star-forming 
dis k of M31 has been completed using the IRAM 30-m 
by iNieten et al.l l)2006l see also references therein). Fi- 
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nally, iMulled (|2003') used the Plateau de Burre interfer- 
ometer to examine the properties of r nolecular clouds in 
9 fields. Using the GAU SSCLUMPS llStutzki fc Giistenl 
ll990tlKramer et al.ll998D algorithm, they decompose the 
emission into 30 individual molecular clouds. 

Previous high-resolution observations of CO in M31 
indicate that a large fraction of the molecular gas 
is found in GMCs. Identifying individual GMCs re- 
quires a telescope beam with a projected size < 5 pc, 
the t ypical size of a CMC in the Milky Way (Blit^ 
^13), which requires an angular resolution of 14" at 
the distance of M31. There have been seven obser- 
vational campaigns that observed CO (1 0) emis- 
sion from M31 at sufficie nt resol u tion t o distinguish 
molec ul ar clouds: Ichikaw a et alJ il985r: IVogel et all 
C987); Lada et al. ( 1988); Wilson & Rudolph" fl993|h 
iLoinard & Allen (1998); Shet h et aL (.2000 ): MuUe] 
With the exception of iLoinard fc AUe n (1998), 
all of these studies have found GMCs with properties 
similar to tho s e fo und in the inner Milky Way and 
'Melch ior et al. | (2000) have argued that the differences 
observed by IToinard & Allen (1998)_ can be attributed 
to observational errors. Indeed. iVogel et al. (1987) pre- 
sented the first direct observations of GMCs in any ex- 
ternal galaxy using interferometric observations. Subse- 
quent studies with interferometers and single-dish tele- 
scopes confirmed that most CO emission in M31 comes 
from GMCs and that the GMCs pro perties were simi- 
lar to those found in the M ilkv Way llLada et alJll98 
IWilson fc RudolDhl l993: Sh eth et al.l2000HMulleil20(r 

Although the molecular gas in M31 has been exten- 
sively studied, there remains a gap connecting the large- 
scale, single-dish observations and the small-scale, inter- 
ferometer observations. To address this gap, we com- 
pleted C0( J =1^0) observations of a large (20 kpc^ 
region) along a spiral arm of M31 with high resolu- 
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tion 50 pc). We then followed up on these obser- 
vations using a more extended configuration of the in- 
terferometer yielding data with a resolution of ~ 25 
pc. This paper presents the observational data of the 
both the survey and the follow-up observations (SjSJ. 
Using only the follow-up data, we present the first re- 
sults, namely a confirmation of previous studies that find 
GMCs in M31 are similar to those in the Milky Way 
(§SJ3EJ- N otably, this paper utilizes the techniques de- 
scribed in ( Ro solowsky fc Ler ov 2006) to correct the ob- 
servational biases that plague extragalactic CO observa- 
tions, thereby placing derived cloud properties on a com- 
mon scale that can be rigorously compared with GMC 
data from other galaxies. The follow-up observations are 
also used to examine the velocity gradients and angular 
momentum of the GMCs, which are then compared to 
the remainder of gas in the galaxy for insight into the 
GMC formation problem We conclude the paper 

by examining the larger galactic environment of M31 to 
explore connections between the GMCs and the larger 
ISM Subsequent work will explore the star forma- 

tion properties of these GMCs and the formation of such 
clouds along the spiral arm using the data from the spiral 
arm survey. 

2. OBSERVATIONS 

We observed ^^CO(J = 1^0) emission from M31 
during the spring and fall observing seasons of 2002 with 
the D and C configurations of t he BIMA millimeter in- 
terferometer fWelc h et al.l ll996^. The observations con- 
sisted of an 81-field mosaic using the most compact (D) 
configuration with follow-up observations on seven sub- 
regions, covering 30 fields at higher resolution (C array). 
The D-array survey spans a projected length of 6.7 kpc 
along a spiral arm in the galaxy. Three of the seven 
follow-up, C-array fields targeted regions with known 
CO emission from the D-array survey, and the remaining 
four fields targeted regions wi th strong CO emis sion in 
the single-dish observations of lDame et al.l l|1993fl over a 
range of galactocentric distances. The locations of the 
fields are indicated in Figure 

The D-array observations were completed in Septem- 
ber and October 2002 over the course of four nights. 
Each night roughly 20 pointings of the mosaic were ob- 
served. During the observations, the fields were observed 
for 46 seconds each, making two passes through the mo- 
saic before returning to the phase calibrator (0102-f 504, 
2.6 Jy) every 30 minutes. This cycle was continued 
through the night, accumulating ~ 6 hours of integration 
time on M31 per night (18 minutes per field). The cor- 
relator was configured to span 200 MHz at 500 kHz (1.3 
km s~^) resolution, easily encompassing the whole range 
of velocitie s expected from t he M31's rotation curve for 
this region l|Ibata et alJl2005l) . At the latitude of the Hat 
Creek Radio Observatory, M31 transits near zenith and 
cannot be observed for 40 minutes a night. This time 
was used to perform flux calibrations using Saturn and 
Uranus. Over the course of the observations, the de- 
rived flux of the phase calibrator was stable to 10%. The 
relative flux calibration of the data may be better than 
this level owing to the intrinsic variability of the phase 
calibrator. 

The data were r educed and inverte d using the MIRIAD 
software package l)Sault et al.lll995|l . following the cali- 
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Fig. 1. — The Digital Sky Survey image of M31 with the loca- 
tions of the interferometer observations indicated. The extent of 
the D-array spiral arm survey is indicated with the black and white 
dashed contour. The targets of the C-array follow-up observations 
are indicated with white regions. Letters next to the C-array fields 
indicate the corresponding fields in Table 111 

bra tion procedures of t he BIMA Survey of Nearby Galax- 
ies llHelfer et al.ll200^ . The four nights of uv data were 
combined using a linear mosaicing technique. The in- 
version used natural weighting, and gridding in the uv 
plane was chosen to produce maps with 3" pixels and 
2.03 km s^^ channel width. The dirty maps were decon- 
volved with a Steer-Dewey-Ito deconvolution algorithm 
optimized for mosaics (M0SSDI2 in the MIRIAD soft- 
ware package). Each plane of the final image was cleaned 
to the l.barms significance level. The final resolution of 
the survey is 14", with negligible variation over the map. 
The mosaic is corrected for the gain of the primary beam. 

Three regions within the D-array survey containing 
emission were observed with the more extended C-array 
configuration. In total, 26 fields were observed on five 
nights from October to December 2002. The C-array ob- 
servations used the same observing strategy as the low- 
resolution survey, observing each field for 46 seconds and 
making two or three passes (depending on the size of the 
mosaic) through the mosaic before observing the phase 
calibrator every 24 minutes. The correlator configura- 
tion and choice of calibrators were the same as for the 
survey. The follow-up observations used a subset of the 
D-array pointing centers to facilitate merging the two 
data sets. These data were combined with the D-array 
observations in the uv plane and inverted using uniform 
weighting. The data were cleaned using the same Steer- 
Dewey-Ito deconvolution algorithm as for the D-array 
map. The final resolution of the combined data is ~ 9" 
with a 2.03 km s~^ channel width. The images are cor- 
rected for the gain of the primary beam. 

Prior to the D-array survey, four fields of M31 were 
observed using the C array (Spring 2002). These obser- 
vations target ed single fields know n to contain bright CO 
emission from lDame et al.l l|1993j) over a range of galacto- 
centric distance. The fields were chosen to span a variety 
of regions in the galaxy (different spiral arms and galac- 
tocentric radii) in an attempt to reveal any biases seen 
by focusing on one region of the galaxy. The observa- 
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tions used 0136+478 as a phase calibrator and Mars to 
establish the flux scale of the phase calibrator. Since the 
observations consisted exclusively of C-array data, they 
were inverted using natural weighting. The final synthe- 
sized beam size for these observations is 7.5". 

The observations produced one survey data cube and 
7 high-resolution data cubes for analysis. The proper- 
ties of these data sets are summarized in Table ^ Since 
space constraints prohibit a full presentation of all the 
observational data, an excerpt from the center of Field 
A is shown in Figure |21 to illustrate typical quality of 
the data. The remainder of the paper emphasizes the 
high-resolution data (Fields A-G in Table P). We defer 
presentation and analysis of the Survey data to a subse- 
quent paper. 

3. ANALYSIS OF THE MOLECULAR EMISSION 
3.1. Signal Identification 

For data cubes A-G in Tabled we identified the emis- 
sion using a dynamic thresholding technique. For every 
pixel in the cube, we first estimated the rms value for 
the noise using the method of Engareiola ct al. ( 2003); 
this method accounts for gain variations in both posi- 
tion and frequency. We first estimate the noise at ev- 
ery position by measuring the standard deviation of the 
data values using iterative rejection of high signal pix- 
els (> 3arms)- We generate a noise map by smoothing 
the resulting pixel estimates with a boxcar kernel of the 
same size as the synthesized beam. Then, we estimate 
the relative changes of the noise across the bandpass by 
estimating arms in each channel map and normalizing 
to the mean value of (Trmsiv) across the bandpass. We 
generate a cube of noise estimates by scaling the error 
map by the fractional change from the mean appropriate 
for that velocity. Performing a pixel-wise division of the 
data cube by the cube of noise estimates yields a data 
set in significance units. 

We identify CO emission in the data cubes by searching 
for contiguous regions of high significance in position- 
position- velocit,j;_S£a£e_jasin£_^ two-tiered thresholding 
scheme (e.g. iScoville et allll987|) . The core pixels in the 
mask are identified as pairs of pixels in adjacent veloc- 
ity channels with I(x,y,v) > 4:armsix,y,v). The mask 
is then expanded in position and velocity space to in- 
clude all pixels with / > 2(7rms that are connected by 
high significance pixels to the Aurms core of the mask. 
This masking process has been used in previous extra- 
galactic, interferometric studies and has been shown to 
succes sfully identify emission with minimal inclusion of 
noise (|Engargiola et al]l2fl03fl . Eliminating noise from 
the emission under consideration is key to the methods 
used to determine cloud properties (see tl3.3|l . 

3.2. Inclination Effects 

The identification and analysis of GMCs in M31 may 
be complicated by the relatively high in clination of the 
galaxy: i = 77° llWalterbos Kenmcuttiil 988^ . The sit- 
uation is similar, in part, to observing GMCs in the outer 
regions of the Milky Way (i = 90°!), and the spatial res- 
olution 25 pc) is comparable to what is found using 
small telescopes (e.g. the CfA 1.2-m) to observe d istant 
(> 10 kpc) GMCs in the Milky Way l|Dame et aU bOOl). 
While study of GMCs has been conducted in the Milky 



Way with such data, the inclination being 77° and not 
90° helps the situation significantly. The foreground and 
background spiral arms in M31 are separated from the 
arm under consideration by > 2' in projection. Thus, 
spatial information is useful for separating the emission 
into GMCs. The clouds in the interarm region are likely 
too faint to be recovered by the deconvolution algorithm 
f ^O) . We conclude that the primary concern in decom- 
positi on is the blending o f GMCs within the same spiral 
arm. Niet en et alJ |2006) estimate that the width of the 
spiral arms in the plane of the galaxy for this region is 
~ 500 pc with a vertical thickness of 150 pc. There could 
be multiple, distinct clouds in the spiral arm along the 
same line of sight, so we make a simple estimate of how 
many GMCs are likely to be along a line of sight through 
the spiral arm. We model the arm as a cylinder with 
elliptical cross-section having major and minor axes of 
500 pc and 150 pc respectively. For a viewing angle 77° 
away from the minor axis of the cylinder, the path length 
through the cylinder is ^ = 400 pc. We assume that the 
GMCs are distributed uniformly in the cylindrical vol- 
ume with volume density n and geometric cross-section 
a. The mean number of clouds intersected by a line-of- 
sight through the arm is then nai. We estimate n by 
counting the number of clouds (N) in some trial volume 
along the arm V = ttLWH/A where L,W and H rep- 
resent the length, width (major axis) and height (minor 
axis) of the volume of the arm respectively. By assum- 
ing a surface density through a typical GMC of Sga/c 
and comparing this to the observed surface density in the 
single-dish map (Eh2)j we can estimate the area filling 
fraction of GMCs in the arm when viewed from above 
the galaxy:/ = (Eh2)/Sga/c- Since / is also given by 
Ncr/LW, we can express 

"'^^^rTT^^ ■ (1) 

Significant confusion will occur for nai > 1. For this 
condition to be met with £ = 400 pc and H = 1 50 pc , 
(I]h2)/Sgji/c > 0.3. The data of iNieten etaP l(200l 
show (Sh2) ~ 5 Mq pc~^ for this region requiring 
GMCs in M31 to have Eqmc < 20 Mq pc'^ for sig- 
nificant blending to occur. We regard this possibility 
as unlikely since GMCs in other, less confused regions 
of M31 are found to have co mparable properties to Lo- 
cal Group clouds llSheth et al . 2000; MuUcr 2003) where 
pnM C « 100 Mq t,c-' ISolomon et a,1in987tiR1itz et alJ 
l2006j) . Unfortunately, one of the aims of the paper is to 
assess how similar clouds in M31 are to those in other 
galaxies; so we cannot reject the possibility that the 
clouds could be significantly blended along the line of 
sight although this would imply significantly larger sizes 
and lower column densities than is typical for GMCs in 
other galaxies. 

We conclude this discussion of inclination effects by 
noting that we will assume that GMC properties are in- 
dependent of viewing angle, as is common in studies of 
GMCs in the Milky Way and beyond (e.g. Scoville ct aH 
^7; Solomon et al. 1987; Wilson & ScoviUc 199^ 
rWilson fc Rudolph J99a: .Sheth et al...200Q) . As a result, 
the reported properties of GMCs do not need any cor- 
rection for the high inclination of the galaxy. 



3.3. Identifying GMCs and Measuring Their Properties 
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TABLE 1 

Summary of BIMA observations of M31. 



r leiQ 




Center 


Resolution 


Size 


Noise 


f 1 

Jrec 






("2000. ^2000) 


(" X km/s) 


(" x" xkm/s) 


(K km/s) 


Survey 


D 


00 41 48.1, +41 19 02 


14 X 2.03 


1885 X 350 X 215 


0.27 


0.35 


A 


C,D 


00 42 26.5, +41 28 45 


8.2 X 2.03 


255 X 205 X 102 


0.64 


0.35 


B 


C,D 


00 42 12.3, +41 25 33 


10.3 X 2.03 


570 X 320 X 102 


0.54 


0.12 


C 


C,D 


00 42 02.4, +41 22 13 


10.2 X 2.03 


305 X 250 X 102 


0.76 


0.17 


D 


C 


00 43 01.8, +41 37 00 


7.0 X 2.03 


100 X 100 X 100 


0.48 


0.42 


E 


C 


00 43 56.8, +41 26 12 


8.1 X 2.03 


100 X 100 X 100 


0.94 


0.23 


F 


C 


00 44 44.3, +41 28 07 


7.2 X 2.03 


100 X 100 X 100 


0.92 


0.15 


G 


C 


00 43 33.9, +41 09 58 


7.1 X 3.04 


100 X 100 X 100 


1.38 


0.19 



^Th e fraction of flux recovered relative to the IRAM 30-m map of the galaxy INieten et alJ 
1200^). 
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Fig. 2. — Channel maps of an excerpt from the center of Field A. Only 12 channel maps containing significant CO emission are shown. 
The gray scale of the channel map runs linearly from to 6 K km s~^. Contours of local significance are shown at —4, —2 (gray), 2, 4, 8, 
and 12 times the local value of arms- At the center of the field, arms = 0.4 K km s~^, though this value changes across the map owing 
to variable gain from the interferometer. Each channel map is annotated in the upper, right-hand corner with Vlsr for the channel in 
km s~^. For reference, the center of the field is at a2000 = 00*^ 42™ 29!3, (S2000 = +41° 28' 53'.'3 and the extent of the displayed region is 



indicated in Figure l3l 

We identify emission in the data cubes from the high- 
resolution data (Fields A-G in Table ^1 using the signal- 
identification method described in Maps of the 
high-resolution data fields appear in Figures 13 and 0| 
This emission was partitioned into individual molecular 
clouds using the segmentation^ algorithm described in 
[Rosolovi'skv & Lerov (2006, RL06) to separate blended 
emission. The segmentation algorithm is a modified wa- 

^ Segmentation is the term used in image processing for dividing 
an image into subregions that share common properties, in this 
case, being part of the same CMC. 



tershed algorithm, related to the original CLUMPFIND 
IWilliams et al. f994). However, the adopted algorithm 
has been optimized to identify GMCs in a wide vari- 
ety of observational data: the algorithm is significantly 
more robust to the presence of noise and does not over- 
segment GMCs as the original CLUMPFIND is known to 
do ( Sheth et al. 2000). The RL06 algorithm estabUshes 
its parameters from physical scales in the data (e.g. par- 
sees, km s""'^) as opposed to observationally determined 
scales (beam and channel widths) to minimize biases in 
comparing catalogs with different observational proper- 
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Fig. 3. — Masked integrated-intensity maps of high resolution fields A-C. The black ellipses indicated the non-deconvolved major and 
minor axes of the molecular clouds identified by the segmentation algorithm. GMCs with sufficient signal to noise to be included in the 
property analysis (Q are shown with a solid black line and low signal-to-noise clouds are shown with dotted ellipses. The gray scale runs 
linearly from K km (indicated with a contour) to the maximum indicated in the upper right-hand corner of each map. Axes are 
aligned with equatorial coordinates and are labeled in offsets measured in arcseconds from the field centers listed in Tables [Tl For reference, 
a scale bar is plotted of projected distance in the plane of the sky and the size of the synthesized beam is depicted beside the field name, 
ff appropriate, the edge of the surveyed region is drawn with a solid line. The square, gray box in the map for Field A indicates the region 
displayed in Figure 1^ 



ties. For the high-resolution data, the algorithm identi- 
fies 67 molecular clouds in the seven cubes. 

For each of these clouds, we measure three macroscopic 
properties: mass, radius (Re), line width (cr^). These 
macroscopic properties are defined using moments of the 
intensity distribution (see RL06 for details). The mass of 
the GMCs is determined by scaling the zeroth moment 
(sum) of the intensity distribution by a constant CO-to- 
H2 conversion factor of 

Xco = 2xl0^° ^ 
K km s 

(jStrong &: Mattoxl Il996t iDame e t al!* '20011 . Assumiirg 

a mean particle mass of 1.36mH (Wilson & Rood 199l| 

gives a mass from the integrated flux over the cloud: 



of the emission distribution along the major and minor 
axes of the cloud: (Jmaj and amin respectively. Like the 
total flux, the values of these moments are extrapolated 
to the K km s~^ pc^ intensity level to minimize the 
effects of the brightness threshold giving cr^^^ and cr^i„ 
(see RL06 for details) . We estimate the deconvolved size 
of the major and minor axes by subtracting the width of 
the beam in quadrature. The size of the cloud is the geo- 
metric mean of the "deconvolved" major- and minor-axis 
sizes: 



beam 



.-,1/2 



2 

beam 



(2) 



Mlum = 4.4 Af( 



o 



CO 



The calculation uses 



K km s 
flux value 



"ipc2' 



W, 



trapolated from the brightness 

using variation m 



CO 



K km s" 



that is ex- 

(2fJr772s) to 

function 



threshold 

[1 Wc.o as 

of th reshold ijScoville et al.l Il987t iRosolowskv fc Lerovl 
200(ill . RL06 demonstrate that this extrapolation prac- 
tically eliminates bias in the data due to low sensitivity, 
though incomplete flux recovery interferometer may still 
affect the resulting mass estimates ( ti3.4(l . An demon- 
stration of the extrapolation method used in this paper 
is shown in Figure [S] 

In a similar fashion to the molecular mass, we use the 
moments of the emission distribution to measure the sizes 
and line widths of the resulting molecular clouds. Fol- 
lowing RL06, the size is defined as the second moment 



The radius of the cloud is estimated by scal ing ctj. 
up by a factor of 1.91 ijSolomon et alJ 119871 S87). 
lBertoldi~fc McKec (199^ note that this radius estimate 
is insensitive to inclination and projection effects. 

The rms velocity dispersion is calculated as the second 
moment of the emission distribution in velocity space. 
Again, the velocity dispersion is extrapolated to to the 
the K km s^^ pc^ intensity level to give the corrected 
velocity dispersion cr*(w). The properties of the molecu- 
lar clouds in the high-resolution data are summarized in 
Table El The clouds are divided into two groups: 'high 
signal-to-noise clouds' with Tpeafc/CTrms > 10 where the 
macroscopic properties can be accurately recovered (see 
> i3.4(l and 'low signal-to-noise clouds' (Tpeafc/CTrms < 10) 
where the derived properties are suspect. 

3.4. Flux Recovery 
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TABLE 2 

Molecular Cloud Properties in M31 



Number Position'' Mlum^ Re^'" Tmax/cTrms I'^v]^ (pv^ 

(",") (10-* Mq) (pc) (kms-l) (kms-lpc-l) (°) 



High Signal-to-Noise Clouds 



1 


(+322, +1246) 


55 


34 


4.07 


17.0 


0.07 


97 


2 


(+340, +1262) 


42 


24 


4.43 


15.1 


0.04 


-151 


3 


(-469, +671) 


72 


52 


6.29 


14.9 


0.12 


-66 


4 


(-238, +828) 


78 


50 


4.86 


14.8 


0.02 


67 


5 


(-300, +764) 


63 


47 


3.38 


14.6 


0.03 


132 


6 


(-256, +811) 


33 


43 


4.44 


13.3 


0.10 


125 


7 


(+2384, +709) 


33 


28 


2.35 


13.2 


0.05 


3 


8 


(-279, +782) 


25 


32 


5.45 


12.6 


0.11 


39 


9 


(+1438, +623) 


31 


19 


2.86 


12.3 


0.03 


106 


10 


(-379, +730) 


34 


33 


4.38 


11.9 


0.06 


43 


11 


(-201, +837) 


15 


24 


2.38 


11.5 


0.06 


58 


12 


(-381, +711) 


61 
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Fig. 5. — Extrapolation of cloud properties (top: non- 
deconvolved cloud radius, bottom: flux) to the i^thresh = 
K km s^'^pc^ flux threshold using the extrapolation algorithm 
of RL06. The extrapolation corrects for signal lost below the 
Scr-ms = 23 K km s~^pc^ clipping level (vertical dotted line). The 
extrapolated value is indicated with a gray dot and the extrapo- 
lation is shown as a gray dashed line. This demonstration of the 
extrapolation is shown for Cloud 19 in Table 1^ which is the lowest 
significance GMC used in the analysis of cloud properties (JD. 

Since interferometers do not measure the total power 
of observed emission, they must rely upon single dish 
observations for accurate measurements of the CO flux. 
In the absence of such single-dish data, deconvolution 
algorithms can be used to extrapolate into the unsample d 
short-spacing area of the (u, v) plane <|Helfer et alJl200^ . 
The accuracy of such extrapolations depends primarily 
on the signal-to-noise ratio of the data and to a lesser 
degree on the observational strategy and the structure 
of the object being observed. Since this work uses flux 
information to calculate cloud masses as well as in the 
moments that determine the cloud properties, flux loss 
affects all stages of this analysis. 

We use the data from IRAM 30-m survey of M31 
(jNieten et al.l |200^ to quantify the amount of flux lost 
in each of the the observations. We report the fraction 
of flux in the BIMA maps relative to that found in the 
IRAM map in Table as the column labeled free- The 
fractions recovered range from 12% to 42%. The rela- 
tively low fraction of flux recovered is a product of two 
factors that affect intcrferome tric observations: spatia l 
filtering of GMCs (studied in ISheth et alllminl l2f)(B 



and the non-linear recovery o f signal in the low signal- 
to-noise reginie_Jemlored inJ^l^^On IIHI^^ The 
work of ISheth et airi)200a i200g) found that, because of 
their clumpy structure, spatial filtering did not hinder 
the accurate measurement of GMC shapes, fluxes and 
properties. However, their simulations were conducted 
in the high si gnal-to-noise regime . In the low signal-to- 
noise regime. iHelfer et al.l l)2002l) found that even rela- 
tively compact simulated sources suffered from signifi- 
cant flux loss, owing to the inability for the deconvolu- 
tion algorithms to isolate low-amplitude signal. These 
two effects (spatial filtering of GMC emission and low 
signal-to-noise) were studied together in RL06 who sim- 
ulated BIMA C,D, and C+B observations of Milky Way 
GMCs (Orion, Rosette, W3/4/5) as if those GMCs were 
located in M31. Their work showed that accurate recov- 
ery of cloud properties required signal-to-noise ratios of 
Tmax/cTrms > 10 for Icss than a 20% loss in all the cloud 
properties. For clouds with lower significance, the flux 
loss is more severe. However, the dominant source of flux 
loss comes from clouds not detected in the observations 
at all, i.e. those with Tmax/ f^rms ^ 1-5, the minimum 
amplitude to which the deconvolution algorithm oper- 
ates. For these clouds, the deconvolution algorithm does 
not identify a significant local maximum for CLEANing 
and the cloud contributes zero flux to the deconvolved 
map. Such low mass clouds likely comprise most of the 
missing flux not seen in the interferometer data. When 
clouds are detected, cloud properties estimated from the 
observations should be relatively accurate, particularly if 
Tmax/cTrms > 10. In general, we only consider the prop- 
erties of clouds in the study if this sensitivity condition 
is met; exceptions to this criterion will be noted. For 
the high-resolution data, 19 GMCs meet the signal-to- 
noise criterion (the 'high signal-to-noise clouds' in Table 
EJ. We emphasize that for these M31 data interferom- 
eter observations act as a filter that preferentially de- 
tects GMCs, but the properties of those clouds are well- 
measured. 

4. MOLECULAR CLOUD PROPERTIES IN M31 

Upon measuring a set of GMC properties for clouds in 
another galaxy, the immediate question is whether these 
clouds are similar to those found in the Milky Way. Ide- 
ally, the well-resolved clouds in the Milky Way would 
be analyzed using the methods of RL06 so that the re- 
sults are directly comparable. However, since measuring 
cloud properties in the inner Galaxy requires establish- 
ing a single distance by resolving the kinematic distance 
ambiguity, a detailed re-analysis has not been completed 
to date. The next best solution is to compare the results 
to an existing catalog. 

For comparison, we use the catalog of iSolomon et al] 
(jl987l S87) to define the population of GMCs in the in- 
ner Milky Way. The work of S87 used similar methods 
to calculate cloud properties though their work adopts a 
fundamentally different procedure for segmenting GMCs. 
In addition, their cloud properties are not corrected for 
the effects of sensitivity as are the properties derived in 
this study. Consequently, there may be systematic dif- 
ferences between their catalog and the present study. In 
general, the fiuxes, line widths and radii of the S87 clouds 
will be smaller than would be derived using the methods 
of RL06; however, using interferometer data to observe 
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the extragalactic clouds will reduce the fluxes, line widths 
and radii by a similar amount (RL06). We have rescaled 
the catalog of S87 to a different galactocentric solar ra- 
dius: Rq = 10 kpc to i?o = 8.5 kpc and the CO-to-H2 
conversion factor adopted in this study. We also recal- 
culated cloud sizes using the definition of RL06. Finally, 
it is noteworthy that the method of RL06 reports, in 
addition to unbiased values for the cloud properties, ap- 
propriate uncertainties for those properties facilitating a 
statistically based comparison of cloud properties in M31 
to those in the Milky Way. 

4.1. Settlings Among GMC Properties 

iLarsonl l|198lD first demonstrated that there were 
power-law scalings among the properties of molecular 
gas clouds. That GMCs followed these scalings became 
well established in a series of surveys and catalogs of the 
CO emitting gas in the inner Milky Way ('Sander s et al.l 
[1985: Da me et a l. 1986; Scovil le et al 1987. S87). Al- 
though there are several algebraically equivalent forms 
for expressing these relations, the most physically rele- 
vant relationships are generally agreed to be: 

and 



For the inner Milky Way, (Jq = 0.6 km s""'^ and a « 1.5 
(S87). The first relationship is the size-line width rela- 
tionship and results from the supersonic turbulent mo- 
tions in the GMCs, as first suggested by Larson (1981). 
The constant a in equa tion 0] is the virial parameter 
l|McKee fc Zweibell IT991 and since a < 2 in the inner 
Milky Way the clouds are frequently interpreted as being 
self-gravitating. In Figure |H1 we plot these two relation- 
ships for the 19 clouds in M31 with sufficient signal to 
noise that their properties can be well-recovered in inter- 
ferometric observations. We compare the relationships 
to those found in the catalog of S87. 

In both the size-line width relationship and the value of 
the virial parameter, there is reasonably good agreement 
between the GMCs in M31 and the inner Milky Way 
GMCs of S87. Fitting the relationship for the 19 clouds 
in M31 that meet the sensitivity criterion f ii3.4|l using 
the method of Akritas & Bcrshadv (1996) gives a size- 
line width relationship of 

logcr„ = (-0.5±0.3) + (0.7±0.2)logi?e. (5) 

For reference, a fit to the catalog of S87 using the same 
fitting method, albeit without uncertainties in the cloud 
properties, gives 

logcr^, = -0.28 + 0.55 logi?e- (6) 

Thus, the size-line width relationship for the molecular 
clouds is indistinguishable from the clouds found in the 
inner Milky Way. The similarity between size-line width 
relationship for clouds in M31 and the inner Milky Way 
has been reported on several occ asions by other authors 
jyoeel et al.' 1987: Wilson & RAidoloh 1993: Gucli n et"an 
l200a .MuUer fc Guehn 2003: Shcth et al. 2000) , though 
no previous study has reported an independent fit for 



the clouds in M31 with well characterized uncertainties 
in the cloud properties. 

The virial parameter for clouds in M31 is constant as 
a function of luminous mass. The uncertainty-weighted, 
mean value of the virial parameter for clouds in M31 is 
(aM3i) = 2.0±0.3. For the sample of S87, (auw) = 1-45. 
The difference between the virial parameters is not sta- 
tistically significant. Moreover, any flux loss will affect 
the luminous mass estimates more than the virial esti- 
mate (RL06), so our measurement of the virial param- 
eter is at worst an upper limit. It is unlikely that the 
true virial parameters of M31 cloud population differs 
significantly from the Milky Way, provided the CO-to- 
II2 conversion factor is the same for both systems. Al- 
ternatively, by assuming that the molecular clouds cata- 
loged in this study are in the same dynamical state as are 
clouds in the inner Milky Way, we can justify our choice 
of CO-to-H2 conversion factor. Again, the similarity of 
the virial paramete r has been noted by several authors 
m previous studies (Lada et al."1988': 'Gue hn et al.1 '2000': 
fuller _& Guehn 2003; Guehn et al. 200|). In particu- 
lar, the latter two studies find no systematic variation 
in the virial parameter for clouds over a large range of 
galactocentric radius in M31. 

5. THE ANGULAR MOMENTUM DEFECT OF GMCS 

5.1. Veloeity Gradients 

In addition to macroscopic properties like line width, 
radius and mass, interferometric observations of ex- 
tragalactic GMCs can be used to measure the ve- 
l ocity gradients a nd angular momenta of the clouds. 
iRosolowskv et alJ (|^D03, RFEB) analyzed the velocity 
gradients of GMCs in M33 to critically evaluate cloud 
formation theories. In this section, we analyze the veloc- 
ity gradients for GMCs in M31 and compare our results 
to those of RFEB and the analysis of the velocity gra- 
dients in GMCs in the inner Milky Way bv iKoda et alJ 

(Unni- 

We determine the velocity gradient by a least-squares 
fit of a plane to the velocity centroid surface Vc{a,S) 
- the velocity centroid measured as a function of po- 
sition across the cloud. The velocity centroid and its 
uncertainty are determined from the first moment of ve- 
locity weighted by brightness. The coefficients of the 
fit determine the magnitude (|Vw|) and position angle 
((/)v ) of the velocity grad ient for the molecular cloud 
(see iGoodman et al.1 119931 for details) . The results of 
the fits are listed in Table El For purposes of this 
analysis, we have included all clouds, since the veloc- 
ity gradient is less affected by the interferometric effects 
which bias the macroscopic properties ft l3.4(l . The veloc- 
ity gradient is derived from line centroids and the first 
moments (centroids) are much less sensitive to the low 
surface brightness wings of the line than are the sec- 
ond moments. We do exclude all measurements of the 
velocity gradient where the uncertainty is larger than 
0.15 km s~^ pc^^. This upper limit was chosen from 
the maximum uncertainty for clouds with high signal-to- 
noise {Tmax/cTrms > 10), including a total of 44 clouds in 
this portion of the analysis. Uncertainties in the derived 
parameters are determined from the errors in the least- 
squares parameters given the uncertainty in the centroid. 

The magnitude of the velocity gradients range between 
and 0.2 km pc~^, in good agreement with both the 
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Fig. 6. — ('e/i) The size-line width relationship for the 19 clouds in M31 with sufficient signal-to-noise for accurate recovery of cloud 
properties. There is good agreement between clouds in M31 and those cataloged in the inner Milky Way by S87. The solid, black and 
dashed, gray lines are fits to the data for GMCs in M31 and the inner Milky Way respectively; the fits are statistically indistinguishable. 
(right) The virial parameter as a function of luminous mass. The virial parameter is roughly constant as a function of luminous mass 
implying a linear scaling between the luminous and virial masses of the clouds. 



values and distributi on found for cfouds in M33 (RPEB) 
and the Milky Way ijKoda et al. 200^. We compare the 
magnitude of the velocity gradient to the velocity gra- 
dient that would be expected from the orbital motion 
of the gas in the galactic potential, ignoring streaming 
motions in the spiral arms. The local galactic velocity 
is estimated from the tangent plane to a surface defined 
by the line of sight projection of the g alactic rotat i on ve- 
locity. We use the rotation curve of llbata et all 
which is a compilation of previous CO and H I studies 
using updated orientation information. On average, the 
rotation gradient of GMCs exceeds the circular galactic 
gradient by a factor of 3.0. 

We measured the difference between the position an- 
gle of the gradient and that of the local galactic gradient 
expected from circular rotation for each of the clouds. 
While suggesting a slight alignment with the galaxy, the 
distribution of this difference is not significantly different 
from a random distribu tion according to a two-sided KS 
test llPress et al.lll992l Pkr = 0.24). This is similar to 
the work of iKodaetahl lj2005r who found a random dis- 
tribution of position angles in the Milky Way GMCs but 
different from RPEB who found significant alignment of 
the gradients with the galactic rotation. 

This analysis has assumed the decomposed data repre- 
sent discrete clouds with independent velocity fields, and 
the factor of three difference between the magnitudes 
of the cloud and local galactic gradients supports that 
model. However, it is also possible that we are inappro- 
priately decomposing a continuous gas flow into clouds 
and the large velocity gradients result from the change 
in the gas flowing into the density wave. In the presence 
of the radial motions that such a flow implies, the local 
galactic velocity gradient will change in magnitude by a 
factor of, at most, \/V^ + Vg /Ve and in position angle 
by ta,n~-^{Vr/Ve) where Vr and Ve are the radial and az- 
imuthal com ponents of th e gas flow. For gas in M31, 
Vr/Ve ~ 0.1 l)Braunlll99H) so these changes will not be 



significant and the factor of three excess of cloud gradient 
magnitude over the local rotational gradient implies dis- 
continuous velocity fields. Moreover, any changes should 
be systematic in nature, and the large range of position 
angles actually observed implies we are measuring the 
velocity gradients of discrete objects. 



5.2. Angular Momentum 

While the magnitudes of the velocity gradients of the 
clouds are larger than the local rotation in the galaxy, 
the implied angular momentum of the GMCs is signifi- 
cantly smaller than that of an equivalent mass of atomic 
gas orbiting in the galactic gravitational potential at the 
GMCs position. The specific angular momentum for 
a cloud of gas with velocity gradient magnitude |Vt;| 
is jGMC = /^IViili^e where P — OA ± 0.1 for a wide 
rang e of density di stributions undergoing uniform rota- 
tion fPhilliDs"'1999') . If the velocity gradient arises solely 
from the turbulen t motions in the cloud, i s likely 2-3 
times smaller ijBurkert fc Bodenheimedl2000() . 

It is unlikely that the high galactic inclination affects 
the results. The near-random distribution of the position 
angles implies that the cloud angular momentum vectors 
are similarly close to random which would require scal- 
ing the angular momentum up by a factor of A/n on 
average (RPEB) . A turbulent origin for the velocity gra- 
dients would not necessitate any correction. Since any 
correction would be small, we make no correction for the 
orientation of the clouds. 

For the clouds in the velocity gradient analysis, the 
mean angular momentum is (36 ± 3) pc km for a 
solid body rotator and 2-3 times smaller for the angular 
momentum of turbulent flows, comparable to the results 
for other systems (RPEB, Koda et al. 206§). We com- 
pare this to the angular momentum of progenitor mate- 
rial for the GMCs. The angular momentum of the gas 
in the galactic potential is given by RPEB (equivalent to 
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where ry is a form factor equal to 1/4 for a contrac- 
tion of a uniform circular disk of material of radius 
ARgai- The length, ARgai, is the range of galactocen- 
tric radii over which material is accumulated to form a 
GMC. The minimum value of ARgai is set by the mass 
of the GMC under consideration and the surface den- 
sity of gas that could comprise the progenitor material: 
Mqmc = T^^ prog^Rlai- Wc take T.prog = 10 M© pc"2 
ijNieten et alJ 12006). including their estimates of both 
atomic and molecular material since small molecular 
clouds seen in CO may be progenitor material for the 
larger GMCs. From the measured masses of GMCs and 
locations of GMCs, we calculate A R„ai and the galac - 
tic shear from the rotation curve of llbata et al.l l)2005f) . 
The mean value of {jgai) = (63 ± 5) pc km s~ for the 
clouds in the sample. On a cloud by cloud basis, the 
angular momentum of the galactic material is only a 
factor of 1.7 larger than the that of the GMCs, signif- 
icantly smaller than the difference found in M33 where 
iigai/jcMc) = 5. The smaller margin is primarily due 
to the lower shear in M31, reducing jgai by a factor of 4 
compared to that found for the clouds in M33 whereas 
the velocity gradients are comparable in the two systems. 
Similarly, the cloud properties and gradients for clouds 
in the Milky Wa y are comparable to the clouds in M31 
ijKoda et al J 1200 5). but the shear is much larger in the 
inner Milky Way compared to the 11-kpc arm of M31. 
Consequently, the discrepancy between the angular mo- 
mentum of progenitor material and the resulting GMCs 
is also larger in the inner Milky Way than it is in M31 (see 
[Blitz 1993). The angular momentum defect for the M31 
clouds is likely significant: the factor of 1.7 difference 
should be regarded as a minimum value since attributing 
the velocity gradients to turbulence or accumulating the 
progenitor material from a larger range of galactocentric 
radius will appreciably increase the margin. 

These results suggest a broad similarity between GMCs 
in M31 and the other disk galaxies in the Local Group. 
The primary goal of studying the velocity structure 
within the GMCs is to relate their angular momentum to 
that of the ISM as a whole to evaluate formation mech- 
anisms. The angular momentum defect between G MCs 
and t he ISM is common among the studied systems filBlitd 
119931 RPEB) and suggests th e importan c e of b raking 
in molecular cloud formation. iKim et all l)2003(l study 
the role of magnetic effects at removing angular momen- 
tum from forming clouds and conclude that magnetic 
effects can reduce the angular momentum to observed 
values. Their unmagnetized simulations show no signifi- 
cant reduction in cloud angular momentum, although the 
density enhancements that represent proto-GMCs show 
jgai/jgmc 2, similar to the observed minimum value 
found in M31. The relatively low-shear environment does 
not necessarily require magnetic braking to produce the 
observed cloud angular momenta, though our measure- 
ment of jgai/jgmc likely significantly underestimates the 
true value. 

6. GALACTIC ENVIRONMENT AND GMC PROPERTIES 



With respect to their macroscopic properties and ve- 
locity gradients, GMCs in the 11-kpc, star-forming arm 
in M31 appears indistinguishable from the molecular 
clouds in the inner Milky Way. As our knowledge of 
GMCs in other galaxies grows, it becomes more apparent 
that such detailed agreement in the properties of GMCs 
is not universal and m ay actually represent exceptional 
cases (iBlitz et al.ll2006|) . Given the differences measured 
among extragalactic cloud populations, the most press- 
ing question is: what physical mechanisms regulate the 
observed differences in GMCs properties? These mecha- 
nisms must couple the GMCs to the ISM on much larger 
scales, and future studies will examine the nature of this 
coupling in more detail. 

The similarity between GMCs in M31 and the Milky 
Way offers a useful point of reference for studying the re- 
lationship between the galactic ISM and GMCs. Finding 
environmental properties that vary significantly between 
the two systems while GMC properties remain the same 
would imply those properties are irrelevant to the regu- 
lation of molecular cloud properties. Several features of 
the galactic environment are thought to play a role in 
regulating the structure of the cold ISM and the proper- 
ties of GMCs. Here, we review several of these features 
and note the degree of similarity to the local Milky Way: 

• Interstellar Radiation Field — The interstellar ra- 
diation field in M31 has been previously noted to 
be poor in the ultravio let relative to the Milky Way 
l|Cesarskv et al.lll998j) . though the M31 field may 
be significantly enhanced in the star-forming spiral 
ar m imaged in the CO observa tions. Indeed, work 
bv iRodriguez-Ferna ndez et al.l (2006) study [C II] 
emission in the star forming arms at Rgai 10 kpc 
and note that the far ultraviolet field is consistent 
with being 100 times the local interstellar value (i.e. 
Go = lO'^) around molecular clouds. Thus, the ob- 
served FUV deficiency likely does not characterize 
the molecule rich region of the galaxy and the radi- 
ation environment is similar to clouds in the Milky 
Way. 

• Metallicity — The gas phase metallicity from H II 
regions in M31 at Rgai = 10 kpc is approximately 
[O/H]-[O/II] (7i=0.0 dex with significant scatter, up 
to ~ 0.2 dex ijDennefeld fc Kunthll98ltlBlair et all 
I1982D . The region shows no signs of being signif- 
icantly enriched or depleted relative to the solar 
neighborhood. 

• Dynamics — The Toomre Q parameter is fre- 
quently used to evaluate the stability of galac- 
tic disks with respect to the formation of self- 
gravitating structures. The parameter is defined 
as 



(8) 



where k is the epicyclic frequency, Cg is the one- 
dimensional gas velocity dispersion and Eg the 
surface density of the gas in the disk. Most 
star -forming regions of galacti c disks show Q « 
1.5 l|Martin fc Kennicuttil200l and this region of 
M31 is no excep tion. For the rotation curve of 
llbata et al.ll|2005fl . k = 0.03 km s"^ pc^^. We take 
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Sg « 10 Mq (|Nieten et al"ll2006j). F inally, we take 
CTg « 8 km from lUnwinI l)1983j) . Combining 
these data give Q w 1.8, well within the range of 
values observed in most galaxies and approximately 
the same value found in the solar neighborhood. 

• Midplane Density — We can also compare the mid- 
plane particle density of this region in M31 to 
the solar neighborhood. The mean particle den- 
sity is important for regulating thermal structure 
of the ISM and the formation of the cold medium 
l|Wolfire et alJl200'^ . Adopting a gas scale height 
of 340 pc for this region of M31 (|Bra,unll199ir) gives 
(n) =0.5 cm~^ in the spiral arm region, but falling 
to (n) = 0.1 cm~'^ outside the arm. These densities 
are smaller than that found in the solar n eighbor- 
hood {{n) = 1.0 cm-3. IWo]fire7nni200l . 

The properties of the ISM in this region of M31 appear, 
at least superficially, similar to those found in the solar 
neighborhood. Hence, the derived similarity of molec- 
ular clouds in the two galaxies is not surprising. The 
only property that appears like it may be significantly 
different between the two environments is the midplane 
particle density in the neutral medium. If the differ- 
ence is significant, the constancy of the GMC proper- 
ties suggests that the particle density may be important 
in regulating the fraction of material found in the cold, 
neutral ISM but perhaps not the properties of individ- 
ual GMCs. While this environmental similarity does not 
isolate any irrelevant physical processes, noting the sim- 
ilarity in both clouds and environments illustrates that 
similar galactic environments do indeed produce similar 
GMCs. 

7. SUMMARY AND CONCLUSIONS 

This paper has presented BIMA millimeter interfer- 
ometer observations of ^^CO(J = 1 — > 0) emission from 
giant molecular clouds (GMCs) along a spiral arm at 
Rgai ~ 10 kpc in M31. The observations were conducted 
in two parts: a low-resolution survey using the com- 
pact (D) configuration of the array and a high-resolution, 
follow-up study using a more extended (C) configura- 
tion of the array. All data from the study were seg- 
mented and analyze d with the analysis algorithms of 
iR.osolowskv fc Lerovl l|20061 . which correct the observa- 
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tional biases imprinted on the data by having relatively 
low signal-to-noise and marginally resolving the GMCs. 
The C-array follow-up spans 7.4 kpc^ finding 67 clouds. 
We report the following conclusions: 

1. The GMC population is very similar to that 
found in the inner Milky Way. Both the size-line 
width relationship and distribution of virial parame- 
ters are statistically indistinguishable from those de- 
rived from the population of clouds in the inner Milky 
Way. Th is confirms the work of several authors 
(iVogel e t al. 1987; Wilso n fc Rudolphll99:jf Guehn et alJ 
l200nt lMullcr & Guclin 2003: She th et alJl200nD while us- 
ing techniques that correct for the biases imprinted by 
using interferometric observations. 

2. We measured the velocity gradients and angular mo- 
menta across the GMCs. Like the Milky Way and M33, 
there is significantly more angular momentum in possible 
progenitor material than in the resulting GMCs, though 
the discrepancy is smallest in M31 owing to the relatively 
low galactic shear at large galactocentric radii. All three 
systems suggest braking mechanisms act during the pro- 
cess of cloud formation to dissipate angular momentum. 

3. The galactic environment in M31 where the GMCs 
are found is similar in most respects that of the solar 
neighborhood and inner Milky Way. That GMCs in both 
environments are similar is consistent with the idea that 
the galactic environment regulates GMC properties. In 
systems where the galactic environment differs signifi- 
cantly from those studied here, the CMC p roperties are 
found to be different re.g.. lBlitz et alJl2006D . 
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